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Abstract 

The  gas  transport  in  the  porous  electrode  is  treated  by  a  phenomenological  approach  such  that  the  gas  concentration  at  the  three-phase 
boundary  (TPB)  region  is  the  additive  superposition  of  that  transported  from  the  source,  i.e.  the  gas  channels.  With  plausible  approximations 
and  elemental  algebra,  analytical  expressions  are  obtained  to  estimate  the  effects  of  ribs  on  the  concentration  polarization  of  planar  fuel  cell 
operations.  It  is  shown  that  the  model  can  closely  reproduce  the  experimental  concentration  polarization  curve  for  small  and  medium  current 
density  (up  to  about  2  A/cm2),  providing  a  simple  and  effective  method  for  engineering  application.  The  concentration  polarization  caused  by 
the  presence  of  a  rib  is  discussed  and  the  concentration  profiles  with  varying  rib  widths  are  illustrated.  In  connection  with  the  electrical 
resistance,  the  determination  of  the  optimal  rib  width  for  minimizing  the  overall  polarization  is  also  shown. 
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1.  Introduction 

Solid  oxide  fuel  cells  offer  great  potential  for  high  effi¬ 
ciency,  environmentally  benign  electric  energy  generation. 
There  is  a  great  deal  of  research  on  solid  oxide  fuel  cells 
generally,  especially  on  new  material  systems  and  on  elec¬ 
trochemical  performance  of  the  cell,  or  membrane-electrode 
assembly  (MEA).  Planar-type  designs  have  received  much 
attention  because  they  potentially  offer  higher  power  density 
relative  to  the  tubular-type  SOFC,  which  is  ascribed  to  the 
low  electrical  resistance  due  to  shorter  current  paths.  In 
particular,  a  very  thin  electrolyte  film  can  be  used  in  an 
anode-supported  SOFC,  drastically  reducing  the  ohmic 
resistance  and  enabling  operation  at  intermediate  tempera¬ 
ture  [1-3].  However,  relatively  little  of  the  literature  con¬ 
centrates  on  such  critical  technical  issues  as  the  system-level 
layout  of  the  layers  and  stacks  [4-6]. 

Under  ideal  conditions  (i.e.  pure  hydrogen  fuel  and  low 
fuel  utilization),  single  cells  with  optimized  materials  can 
yield  power  densities  on  the  order  of  0.5-1  W/cnr  or  higher. 
For  practical  applications,  however,  multiple  cells  must  be 
connected  in  series  by  interconnect  material  to  form  a  fuel 
cell  stack.  Small  channels,  which  are  often  formed  in  the 
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interconnect  material,  are  commonly  used  to  carry  the  fuel 
and  air  gas  flow.  The  “ribs”,  which  separate  and  define  the 
flow  channels,  make  direct  contact  with  the  electrodes.  In 
designing  the  layer  architecture,  there  is  a  tradeoff  that  must 
be  considered  between  the  rib  and  channel  sizes.  On  one 
hand,  wider  ribs  and  ribs  covering  bigger  fraction  of  the  cell 
area  may  reduce  the  interface  resistance  to  current  flow  by 
increasing  the  electrode-interconnect  contact  area  and  redu¬ 
cing  the  current  path  through  the  possibly  high  resistance 
electrode  material.  Hence,  such  ribs  will  give  a  better 
conduction  of  the  electrical  current  and  reduce  ohmic  losses. 
On  the  other  hand,  the  chemical  species  do  not  diffuse  as 
well  underneath  wide  ribs.  Narrow  ribs  are  needed  to  facil¬ 
itate  more  uniform  distribution  of  reactive  gases  across  the 
area  of  the  electrolyte  surface  and  thus  to  promote  electro¬ 
chemical  performance.  The  implication  of  the  tradeoff  to 
the  cell  performance  is  very  significant.  It  was  reported  that 
the  maximum  power  output  was  only  87 1  mW/cm2  for  the 
combination  of  4  mm  channels  and  4  mm  ribs,  while  that  for 
1  mm  channels  and  2  mm  ribs  was  1930  mW/cnr  for  some 
proton  exchange  membrane  (PEM)  fuel  cells  [7].  Clearly,  a 
designer  must  decide  how  to  configure  a  channel  network  to 
optimize  the  cell  performance. 

The  general  approach  in  discussing  the  rib  effects  is  to 
map  out  the  transportation  process  in  the  case  of  rib  pre¬ 
sence,  starting  from  the  known  parameters  when  ribs  are 
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absent.  A  rigorous  analysis  requires  a  3D  model  of  the 
porous  electrode  joined  with  the  interconnect  rib  and  elec¬ 
trolyte  sheet.  The  gas  transportation  path  determines  the 
diffusion  polarization  and  electrochemical  performance,  and 
the  electrical  path  determines  the  ohmic  polarization.  The 
optimal  rib  design  can  then  be  determined.  While  such  3D 
analysis  is  undoubtedly  useful  in  revealing  the  details  of  the 
complex  interaction  of  competing  factors,  the  results  are 
difficult  to  generalize  because  they  are  usually  dependent  on 
a  great  number  of  design  and  operational  parameters.  As 
such,  the  approach  is  usually  unable  to  provide  easy-to-use 
design  guidance  to  the  broad  fuel  cell  engineering  society. 

The  purpose  of  the  present  work  is  to  derive  some  simple 
guidelines  for  estimating  the  effects  of  the  interconnect  rib 
area  on  the  concentration  polarization  as  a  function  of  the 
different  physical  and  geometrical  parameters  involved. 
This  will  provide  a  more  complete  basis  for  optimizing 
the  stack  repeat  unit  geometry — especially  the  shape  and 
dimensions  of  the  interconnect.  The  results  of  this  paper 
assist  the  designer  in  understanding  how  channel  sizing 
affects  gas  uniformity  for  the  system  by  formulating  the 
model  in  a  dimensionless  framework,  enabling  a  general¬ 
ization  of  the  results  in  terms  of  one  characteristic  diffusion 
parameter  and  one  dimensionless  geometry  parameter. 

2.  Concentration  polarization  model 

The  maximum  usable  work  at  constant  temperature 
and  pressure  is  the  free  energy  change  of  a  reaction, 
G  =  Go  +  ( RT /2F)  \n(Pu2Plf  / Pu2o),  where  G0  is  the  free 
energy  change  at  the  standard  state  (for  simplicity,  we  use  H2 
fuel  in  the  discussion,  though  the  principle  is  of  general 
applicability).  It  can  be  easily  seen  that  the  fuel  gas  pressures 
contribute  to  the  Nernst  potential  by  [RT /2F]  I n |/>[ |2 //jh2cj1. 
At  gas  channels,  the  fuel  gas  pressures  are,  respectively,  P^ 
and  o-  These  correspond  to  the  maximum  work  available 
for  external  work  if  the  electrochemical  reaction  were  taking 
place  at  the  channels;  however,  the  reaction  takes  place  at  the 
three-phase  boundary  (TPB)  (i.e.  gas-electrolyte-electrode). 
The  gas  pressures  at  TPB  are,  respectively,  P\i,  and  Ph2o  due 
to  diffusion  resistance  (loss).  The  maximum  work  available  is 
then  only  Go  +  RT\n\P-a2/Pn,c)[,  and  the  Nernst  potential  is 
less  than  the  ideal  case  by  [RT /2F]  hi{PH,PH2o/f3H2f3H2o}- 
The  loss  is  due  to  the  anode  concentration  polarization. 
Similarly,  resistance  to  the  oxidant  diffusion  from  cathode 
gas  channel  to  the  corresponding  TPB  region  causes  the 
cathode  concentration  polarization. 

Suppose  a  point  X  at  the  TPB  has  coordinate  ( x ,  y,  0),  a 
point  Q  is  at  gas  channel-electrode  interface  with  (m,  v,  Ze), 
where  /e  is  the  electrode  thickness.  Without  loss  of  general¬ 
ity,  we  can  assume  the  channels  are  along  the  y-direction.  Let 
L0  be  the  penetration  distance,  characterizing  the  maximum 
distance  the  gas  can  diffuse  from  the  channel  before  it  is 
consumed  by  the  cell  reaction.  From  a  phenomenological 
point  of  view,  the  gas  concentration  at  X  is  the  summation  of 


those  diffused  from  sources  of  all  <2’s  at  the  electrode¬ 
channel  surfaces  and  within  the  penetration  distance  from  X: 

Px  =  JdPx  =  J  C(u,v)dAQ,  (1) 

where  cL4q  =  d u  dv  is  the  area  represented  by  Q,  and  the  3D 
diffusion  source  function,  C,  is  expected  to  be  dependent  on 
the  gas  pressure  at  Q  and  the  distance  between  Q  and  X. 
Integration  along  the  y-direction  (channel)  can  be  carried  out 
implicitly  by  using  average  property  of  C  and  reducing  the 
3D  diffusion  form  to  an  in-plane  diffusion  expression: 

/rx+Lo 

d Px=  /  B{u)F{u)du  (2) 

Jx-Lo 

Here,  B(u)  is  the  in-plane  source  function,  while  /*’(;<)  takes  a 
value  of  1  when  u  is  inside  channels  and  0  when  a  is  inside  the 
ribs.  It  can  be  expected  that  B(u)  decays  roughly  exponen¬ 
tially  with  the  distance  between  X  and  Q  due  to  the  nature  of 
concentration  diffusion.  One  can  quite  reasonably  write 

B  =  at~cdi  (3) 

so  that  the  gradient  of  the  source  function  is  balanced  by 
the  flow  resistance  (proportional  to  the  current  density)  in 
accordance  with  Darcy’s  law  of  diffusion  in  a  porous 
medium  [8].  Here,  c  is  a  constant  characterized  by  the 
electrode  diffusivity,  best  determined  by  the  experimental 

data  and  cl  =  'll  +  {u  —  x)~  is  the  distance  between  u  and 
x.  The  effective  current  density  in  the  channel  is  i,  and 
i  =  *o/ ( 1  — /),  where  i0  is  the  current  density  without  a  rib 
present,  and /is  the  fraction  of  rib  width  to  the  rib-channel 
width,  as  seen  by  the  point  X.  That  is,  the  effective  channel 
gas  flux  (proportional  to  i)  should  be  1/(1  — /)  times  that  of 
no-rib  flux  due  to  smaller  cross-section  available  for  gas 
transport  with  rib  presence.  Assuming  no  capillary  forces, 
Px  =  Po,  when  i  =  0,  and  it  is  found  a  =  Pq/ [2Lq  (1  —  /)] .  In 
other  words,  the  source  function  is 

,4> 

By  integrating  Eq.  (2)  for  the  anode,  one  can  obtain  the  H2 
pressures  at  the  TPB  with  a  single  undetermined  parameter, 
“c”,  when  the  electrode  thickness  and  the  current  density 
are  specified.  The  best  fit  to  the  experimental  data  without 
rib  presence  provides  the  desired  “c”.  With  the  specified 
“c”,  integration  of  Eq.  (2)  with  rib  presence  would  allow  us 
to  discuss  the  rib  effects  on  concentration  polarization. 

Applying  a  similar  approach  to  the  cathode  side,  one  can 
obtain  02  pressure  at  the  TPB.  Notice  that  there  is  in  general 
no  correlation  between  “c”  for  H2  and  “c”  for  02  due  to 
different  microstructures  and  materials.  However,  it  is  rea¬ 
sonable  to  expect  “c”  to  be  inversely  proportional  to  the 
diffusion  coefficient,  and  one  may  set  c(02)  =  (1/2) 
c(H2)[D(H2-H20)/D(02-N2)][e(H2)T(02)/£(02)T(H2]« 
2[£(H2)t(02)/£(02)t(H2]c(//2).  The  first  factor,  1/2,  is  due 
to  only  half  the  amount  of  02  as  of  H2  needing  to  be 
transported.  Here,  e/t  is  the  porosity/tortuosity  factor;  their 
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inclusion  changes  the  factor  “2”  to  some  different  value  and 
does  not  change  the  qualitative  results  described  below  for 
anode-supported  fuel  cells,  where  cathode  diffusion  polar¬ 
ization  is  small  compared  with  that  of  the  anode  and  will  not 
be  discussed  explicitly. 

Eq.  (2)  generally  requires  numerical  integration.  Though 
numerical  integration  poses  no  difficulties  in  principle,  an 
analytical  approach  is  preferred  for  ease  of  analysis  and 
is  discussed  below.  Assuming  the  electrode  thickness  is 
relatively  small  compared  with  L0,  we  may  approximate 
d  =  [I2  +  (u  —  x)2]  ~  \u  —  x\.  This  approximation  is  reason¬ 
able  for  a  thin  electrode  or  low  to  medium  current  density, 
where  the  penetration  distance  is  expected  to  be  large,  as  it  is 
in  the  case  we  are  focused  on.  For  high  current  density,  such 
approximation  is  not  expected  to  be  very  good;  however,  we 
will  show  that  this  is  not  a  serious  limitation  for  practical 
purposes. 

When  there  is  no  rib  present,  one  can  integrate  Eq.  (2) 
easily  if  d  is  replaced  with  \u  —  x|  and  gives  the  TPB  H2 
pressure  as 


/y  =  p 


H  2‘ 


(1  -  e~A°) 


(5) 


where  Ao  =  cLoio. 

For  H20  pressure  at  the  TPB,  which  is  higher  than  that  at 
the  channel-electrode  surface  for  H20  transport  out  of  the 
TPB,  a  treatment  similar  to  that  for  H2  can  be  used.  The 
resulting  expression  depends  on  Ch2o>  which  is  in  principle 
different  from  that  for  H2  due  to  different  diffusivities. 
Because  H2  is  lighter  and  moves  faster  than  H20,  the  total 
pressure  of  H20  at  TPB  should  be  slightly  higher  than  that  at 
the  channel-electrode  surface  for  an  equal  counter-current 
of  H2  and  H20  to  occur.  However,  this  adds  undue  complex¬ 
ity  with  an  additional  parameter  and  is  not  necessary  for  our 
discussion,  which  is  to  provide  easy-to-use  guidelines  for 
estimating  the  rib  effects.  Therefore,  we  will  simply  follow 
Kim  et  al.  [2]  and  assume  P\\2  +  Ph2o  =  P,  where  P  is  the 
total  pressure  of  H2  and  H20  in  the  gas  channel.  Therefore, 
the  TPB  H20  pressure  can  be  obtained  by 


h2o 


—  p° 

—  rH20 

^  /j0h2  -  /y2 

II 

*0 

K° 

O 

1+  P°H 2 
p° 

[  ru2o  L 

1  - 


(1  -  e~A°) 


(6) 


The  anode  concentration  polarization  is  then 
//(diffusion,  Ao) 

(1  -  e~A°)/Ao 


1  +  (^h2/^h2o)[1  -  (1  -  e-^o)/A0] 


(7) 


One  can  find  the  parameter  c  by  fitting  Eq.  (7)  to  experi¬ 
mental  anode  polarization  data.  Rewriting  the  usual  anode 
polarization  form, 


/  RT  I 

^(diffusion)  =  -  —  In 


1  (*/4s) 


1  + W^oXi/ias), 


(8) 


2 


*2 


cs 

'o 

pH 


rs 

© 

< 


Current  Density(A/cm2) 


Fig.  1.  Comparison  of  present  concentration  polarization  model  with 
Eq.  (8)  and  experimental  parameters. 


Virkar  et  al.  [1]  and  Kim  et  al.  [2]  found  the  best  fit  to  their 
experimental  data  as  i-ds  =  4.74-6.35  A/cm2  for  intermedi¬ 
ate  working  temperature  ( T  =  700—800  °C),  where  zas  is  the 
anode  limiting  current.  Although  the  anode  polarization 
expression  of  Eq.  (7)  is  very  different  from  that  of  Eq.  (8), 
using  A// =  0.36  and  0.51  (cm2/ A),  respectively,  for 
4s  =  4.74  and  6.35  A/cm2,  Eq.  (7)  can  fit  Eq.  (8)  reasonably 
well.  Fig.  1  shows  the  comparison  of  Eqs.  (7)  and  (8)  for  the 
two  representative  cases.  As  can  be  seen,  the  fit  is  good  for 
small  and  medium  current  density.  The  difference  is  less 
than  5  mV  for  current  density  up  to  about  2  A/cm2.  For 
larger  current  density,  the  fit  is  not  good  due  to  a  much  faster 
increase  of  Eq.  (8)  than  that  of  Eq.  (7).  This  is  a  consequence 
of  the  approximation  made  in  calculating  distance,  d,  where 
the  anode  thickness  is  neglected.  However,  it  is  not  a 
serious  limitation  because  the  fit  is  good  for  working  current 
densities.  Therefore,  Eq.  (7)  is  a  useful  expression  in 
describing  anode  polarization. 


3.  Rib  effects 

Due  to  the  smaller  cross-section  available  for  gas  trans¬ 
port,  the  ribs  lengthen  the  effective  gas  diffusion  distance, 
resulting  in  increased  concentration  polarization.  In  parti¬ 
cular,  if  the  rib  width  is  larger  than  the  diffusion  penetration 
distance,  there  will  be  appreciable  underused  or  even  unused 
areas  under  the  rib  shadows,  reducing  the  cell  performance 
substantially.  On  the  other  hand,  if  the  rib  width  is  smaller 
than  the  penetration  distance,  the  gas  concentration  at  the 
TPB  is  expected  to  be  uniform,  and  the  cell  should  perform 
relatively  well.  Therefore,  from  the  concentration  polari¬ 
zation  point  of  view,  narrow  ribs  are  preferred  over  wide 
ribs  if  the  ratio  of  rib  width  to  that  of  channel  width  is 
fixed. 
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3.1.  Narrow  rib  effect 


Assuming  uniform  channel  spacing,  the  rib  width  is  a,  and 
the  channel  width  is  cl  —  a;f  =  aid  is  the  fraction  of  the  rib 
width  in  the  width  of  rib-channel  pair,  d.  Without  loss  of 
generality,  we  may  choose  the  TPB  location  as  0  <  x  <  d.  If 
x  <  a,  the  point  is  under  the  rib  shadow,  while  that  for  x  >  a 
is  under  the  channel.  The  left-most  channel  within  the 
penetration  distance  is  at  N \  <  0;  the  right-most  channel 
is  at  N2  >  0  ( N2  is  expected  to  be  |Afi|  or  |/V|  |  ±  1).  The 
channels’  ranges  are  [jd  +  a,  (j  +  l)d],  and  the  ribs’  ranges 
are  (JdJd  +  a),j  £  [TV] ,  Nf\.  For  cases  with  very  narrow  ribs, 
Lo  =  Nd  =  Nid  =  —N\d  S>  d,  Eq.  (2)  can  be  written  as 

rd 

s-A(\jd+u-x\/Lo) 


i\  = 


Po 

.2X<,(1  -/) 


N2  ra 

w 


j=Ni“ 


2Lq  ( 1  -/) 


j=-NJ 


QA(jd+u-x)/Lo 


+ 


N 

E 


-A(jd+u—x)/Lo  du  . 


-A\u-x\/L0  du 


(9) 


After  some  algebra,  we  have 


Pr  = 


Po 

2A(l  —  f) 

~(x/d)(A/N) 


x  e 


1  -  &~AIN 

+  e-((d-x+a)/d)(A/N^  +  ^ 


U-((.“-x)/d)(A/N)  _  e- 


((d-x)/d)(A/N)}  +  §(x  _ 


x  [2  —  e 


(( a-x)/d)(A/N ) 


—  e 


-((d-x)/d)(A/Nh 


(10) 


Considering  very  narrow  rib  limit  (N  8>  1),  Eq.  (10)  reduces 
to 

(1-0 


Ph2  =  P 


h2 


(ID 


Fig.  2.  Rib  contribution  to  the  concentration  polarization  as  a  function  of 
the  rib  width  fraction,  /. 

Fig.  2  shows  the  extra  concentration  polarization  with 
/  =  0.2,  0.4,  0.6,  and  0.8  for  various  current  densities. 

For  anode-supported  fuel  cells,  the  cathode  concentration 
polarization  is  expected  to  be  much  smaller  than  that  of  the 
anode.  The  extra  concentration  polarization  caused  by  ribs  in 
the  cathode  side  is  also  expected  to  be  small  and  is  not 
discussed  explicitly  here. 

3.2.  Wide  rib  effect 

When  the  penetration  distance  is  about  half  the  rib  width 
or  smaller,  the  rib  effect  on  the  concentration  polarization  is 
expected  to  be  the  most  significant  due  to  significant  under¬ 
utilized  area  under  ribs.  In  the  following,  we  focus  on  cases 
with  Lq  <  d.  In  such  cases,  the  possible  integrals  involved  in 
calculating  Eq.  (2)  are 


where  A  =  cL0i  =  cL0i0/(l  -/)  =  A0/(l  -/).  This 
expression  has  the  same  form  as  the  no-rib  case  except  that 
the  effective  current  density  should  be  used.  In  other  words, 
the  no-rib  expression  is  a  special  case  of  the  narrow  rib 
expression  with  /  =  0.  The  gas  concentration  is  uniform  in 
the  TPB  region  whether  the  area  is  under  channels  or  under 
ribs.  The  uniformity  is  due  to  the  diffusion  penetration  dis¬ 
tance  being  much  larger  than  the  rib  width.  In  fact,  numerical 
results  show  that,  in  practice,  uniformity  is  achieved  with  any 
N  >  1.  Therefore,  Eq.  (11)  is  a  representative  expression  for 
all  cases  where  the  diffusion  penetration  distance  is  larger  than 
the  combined  width  of  a  rib  and  a  channel. 

As  the  above  derivation  shows,  the  net  effect  of  narrow 
ribs  is  to  replace  A0  with  A  =  Ao/(  1  —  /)  in  Eq.  (5),  where/ 
is  the  fraction  of  rib  width  to  the  width  of  the  rib-channel 
pair.  The  difference  between  using  A  and  A0  in  Eq.  (7)  is  the 
extra  anode  polarization  caused  by  the  narrow  ribs: 

J7rib  (diffusion,  A)  =  (/(diffusion,  A)  —  (/(diffusion,  Ao) 

(12) 


r  0 

h=  qA(ii-x)/Lo  du  §(1^  _  x)S(x  —  Lq  +  d  —  a) 

J  x—Lo 

=  —  (e~A'/i°  —  e~A)(5(Lo  —  x)<5(.r  —  Lq  +  d  —  a)  (13a) 
A 

Here,  the  step  functions  are  used  to  indicate  that  /,  appears 
only  if  x  <  Lq  and  x  —  Lq  >  —d  +  a. 

ro 

I2=  e-A\u-x\/Lo  dM  §(_d  a  —  X+  Lq) 

J  —d+a 

=  —  [e-^  -  ^d-a+xP^]5{-d  +  a-x  +  Lq)  (13b) 
A 

rd 

h=  erA\u~x\/L<'  d«  5(a  —  x  +  Lq)5(x  +  Lq  —  d) 

J  a 

=  —  {6{a  -  x)  [&-A(p-*)/Lo  -  e-A(d-x)/L 0] 

A 

+  5(x  -a)[ 2  -  eA(fl-*)/io  - 
x  5(a  —  x  +  Lq)S(x  +  Lq  —  d) 


(13c) 
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rx-j-Lo 

h  =  e-A\u-x\/L0  dM  _  x  _|_  Lq)5(x  +  Lq  —  a) 

J  a 

xd(d-x-Lo)  =  {<5(«  -  x) -  e~A] 

+  S(x  -a)[ 2  -  eA(a-*)/I°  -  e-A]}<5(a  -  x  +  L„) 
x  8(x  +  Lq  —  a)d(d  —  x  —  Lq)  (13d) 

rd 

h=  erA\u~x\lL(>  &u  8(x  —  Lq  —  ci)d(x  +  Lq  —  d) 

Jx-Lo 

=  —  (2  -  e^”*^  -  e~A)(5(x  -  L0  -  a)<5(x  +  L„  -  d) 

A 

(13e) 

rx+Lo 

h=  e~A\u~x^Lo  du  S(x  —  Lo  —  a)S(d  —  x  —  Lo) 

Jx-Lo 

=  —  (2  -  2e~A)<5(x  -  Lo  -  a)<5(d  -  x  -  L0)  (13f) 

A 

/>*+Lo 

Ij  =  I  e-A|«-jc|/io  <5(x  +  L0  —  d  —  a) 

J  d+a 

=  ^  ^-A(d+a-x)IU,  _  e-A)^(x  +  L0  —  J  —  a)  (13g) 

/l 

With  the  above  integrals,  the  gas  pressure  is  found  by  simple 
summation: 

7 

Px  =  yYjj 

7=1 

Despite  their  lengthy  appearance,  the  above  analytical 
expressions  contain  only  elemental  functions  and  can  be 
easily  implemented  in  a  spreadsheet.  As  an  example.  Fig.  3 
shows  the  concentration  distribution  for  Lq  =  d/A  and  the 
rib  width  =  (0.2,0. 4, 0.6)<iforthe  current  density  of  1  A/cm2. 
The  concentration  distributions  with  no  ribs  and  with 
narrow  ribs  are  also  shown  for  comparison.  The  average 
concentrations  for  the  rib  widths  of  0.2,  0.4,  and  0.6  are  98, 


Fig.  3.  Fuel  concentration  distribution  for  Lo/d  =  1/4  with  various  rib 
widths. 


87  and  74%,  respectively,  of  that  for  the  corresponding 
narrow  ribs.  Clearly,  wide  ribs  cause  higher  concentration 
polarizations,  and  designs  with  narrow  ribs  are,  in  general, 
preferred  if  the  percentage  of  rib  width  is  fixed. 

4.  Rib  optimization 

The  optimal  rib  designs  depend  on  balancing  the  rib 
effects  on  the  concentration  polarization  and  the  electrical 
resistance.  The  electrical  resistance  includes  both  the  ohmic 
resistance  of  the  rib  material  and  the  contact  resistance 
between  the  ribs  and  the  electrodes.  Typically,  the  contact 
resistance  is  much  higher  than  the  ohmic  resistance  of  the  rib 
material  [7,9-13].  The  overall  area  electrical  resistance  of 
the  ribs  can  be  written  as  R =  /(//,  where/is  the  rib  width 
fraction  as  discussed  above  and  typical  /(-values  are  10- 
100  mQ  cm2  [7].  Without  the  contact  resistance,  [>  is  usually 
smaller  than  1  mQ  cm2  [9,10].  For  simplicity,  we  will  only 
discuss  the  narrow  rib  scenario.  More  general  discussion  and 
characterization  of  the  rib  design  will  be  reported  in  another 
paper. 

Fig.  4  shows  the  variation  of  the  overall  overpotential 
caused  by  ribs  with /for  three  representative  /-values  of  2, 
10,  and  50  mQ  cm2  for  current  density  of  1  A/cm2  and 
operating  temperatures  of  800  °C  and  /as  =  5.65  A/cm2. 
The  minimum  overall  polarization  indicates  the  best  overall 
rib  width  design  for  the  specified  operating  temperature  and 
target  current  density.  Though  small  rib  width  fraction  is 
intuitively  favored,  the  optimal  rib  width  depends  strongly 
on  the  electrical  resistance.  The  optimal  rib  width  fraction  is, 
respectively,  /  =  0.2,  0.35,  0.6,  for  /  =  2,  10,  50  mQ  cm2. 
The  optimal  rib  width  is  similar  for  T  =  800  °C  and 
T  =  700  °C.  The  higher  the  electrical  resistance  caused 
by  ribs,  the  larger  is  the  rib  width  fraction  preferred.  Small 
rib  width  fraction  (f  <  0.2)  is  preferred  only  when  the 
contact  resistance  is  very  small.  In  most  cases,  the  rib  width 
fraction  of  >1/3  is  expected. 


Fig.  4.  Overall  rib  polarizations  by  varying  rib  width  fractions  for  three  rib 
electrical  resistances. 
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5.  Summary 

We  have  provided  a  simplified  phenomenological  model 
for  the  concentration  polarization.  The  model  fits  the  experi¬ 
mental  data  well  for  any  reasonable  working  current  density. 
Using  the  model,  the  effect  of  rib  size  on  the  concentration 
polarization  can  be  obtained  analytically,  and  the  model 
provides  an  easy-to-use  guidance  for  optimizing  the  rib- 
channel  layout.  The  following  conclusions  are  based  on  the 
results  from  the  analytical  solutions: 

•  When  the  rib  width  is  small  compared  with  that  of  the 
characteristic  penetration  distance,  the  gas  concentration 
is  uniform.  The  concentration  polarization  with  the  rib 
presence  takes  the  same  form  as  that  without  ribs  except 
that  the  effective  current  density  should  be  used. 

•  When  the  rib  width  is  comparable  or  larger  than  that  of  the 
characteristic  penetration  distance,  non-uniform  gas  con¬ 
centration  results  and  the  average  concentration  is  smaller 
than  that  of  the  narrow  ribs  for  the  corresponding  rib 
width  fraction,  resulting  in  higher  concentration  polariza¬ 
tion.  Therefore,  narrow  ribs  are  preferred  over  wide  ribs 
for  a  chosen  rib  width  fraction. 

The  optimal  rib  design  is  obtained  by  minimizing  the 
overall  ohmic  and  concentration  polarization  of  the  ribs.  The 
optimal  rib  width  is  affected  significantly  by  the  electrical 
resistance.  The  higher  the  rib  electrical  resistance,  the  larger 
the  rib  width  faction  should  be.  For  realistic  electrical 
resistance,  the  rib  width  fraction  is  expected  to  be  between 
1/3  and  2/3  of  the  channel  width. 
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